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[ Abstract. We derive a uniform exponential decay of the total energy for the 

^ ' nonlinear Klein-Gordon equation with a damping around spatial infinity in or 

in the exterior of a star-shaped obstacle. Such a result was first proved by Zuazua 
[401 141j for defocusing nonlinearity with moderate growth, and later extended to 
the energy subcritical case by Dehman-Lebeau-Zuazua 7 , using linear approxima- 
P I ' tion and unique continuation arguments. We propose a different approach based 

I solely on Morawetz-type a priori estimates, which applies to defocusing nonlinear- 

ity of arbitrary growth, including the energy critical case, the supercritical case 
and exponential nonlinearities in any dimensions. One advantage of our proof, 
even in the case of moderate growth, is that the decay rate is independent of the 
nonlinearity. We can also treat the focusing case for those solutions with energy 
less than the one of the ground state, once we get control of the nonlinear part 
in Morawetz-type estimates. In particular this can be achieved when we have the 
^ ' scattering for the undamped equation. 
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2 l. aloui, s. ibrahim, and k. nakanishi 

1. Introduction 

We study total energy decay for the damped nonlinear Klein-Gordon equation 
(NLKG) for u{t, x) : [0, oo) x ^ R, 

Utt + a{x)ut — Au + U + f'{u) = (t, x) G [0, oo) x Q, 

u = {t,x) e {0,oo) X dn, ^^'^"^ 

where Q is a (exterior) domain in (A^ > 1), the damper a : f2 — > [0, oo) is 
effective uniformly around the spatial infinity 

a{x)>0, a G L°°(fi), liminf a(x) > 0, (12) 
and the nonlinear energy / : M — )■ M is satisfying 

f{u) = f{\u\), m = /'(o) = no) = 0. (1.3) 

This equation may be regarded as a simple model for nonlinear dissipative hyperbolic 
equations with purely dispersive spatial regions. The stabilization problem consists 
in deriving sufficient conditions on the dissipation a{x) guaranteeing uniform decay 
of the total energy, defined by 

E{u;t) := / \ut\'^ + \Vu\'^ + \u\'^ + 2f{u)dx. (1.4) 
Jq 

When the nonlinear energy is not positive definite in the whole energy space, we 
have to distinguish it from the free energy 

EF{u]t):= / \u\^ + \Vu\^ + \u\^dx. (1.5) 

Multiplying (11. ip by ii and integrating over the slab (^1,^2) x one formally gets 
the energy identity 

E{u;t2) - E{u;ti) = -2 [ [ a{x)\u{t,x)\'^dxdt. (1.6) 



Heuristically, one achieves the stabilization when certain portion of the energy keeps 
flowing into the kinetic part on the region of effective dissipation. Then one typically 
obtains exponential decay in the form 

E{u;t) <Ce-^'E{u;0), (1.7) 

with some constants C, 7 > 0. 

The two major factors in this problem are the geometric structure of the domain, 
the boundary condition and the dissipation on one hand, and the nonlinear inter- 
action on the other hand. The geometric aspects have been extensively studied 
especially in the linear setting, with deep insights into wave propagation properties, 
see e.g., [BEHIE]. 

However the nonlinear interaction becomes the bigger obstruction for the finer 
investigation of the wave propagation. The standard strategy is to start with some 



""^These assumptions are quite natural for the nonlinear Klein-Gordon energy. /(O) = is 
necessary to have finite energy, /'(O) = follows from the evenness, while /"(O) = means that 
the mass term is exactly u. It can be replaced with mu (m > 0), but in the massless case m = 
we do not have the exponential decay. See e.g. [27] for the algebraic decay in the massless case. 
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fairly weak decay properties which are accessible even in the nonlinear setting, and 
then through some limiting procedures and approximation arguments, to reduce the 
problem to the linear case, where much finer informations are available. In this type 
of arguments, the linear approximation poses essential limitations on the generality 
of the nonlinear interaction. We refer to |10l SH |5l [TJ |6l [1] for the stabilization 
results in the current setting (11. ip . 

1.1. Defocusing case. Our observation is that one can skip the latter part of the 
argument, i.e. the linear approximation, at least in the simplest geometric setting, 
thereby obtaining the exponential decay for fairly general nonlinearity with the 
repulsive sign. Our nonlinearity condition is a sort of self-coerciveness, given in 
terms of 

g{u):=uf'{u)-2f{u). (1.8) 
The first main result of this paper is 

Theorem 1.1 (defocusing case). Let Q C (N & N) be either the whole space or 
the exterior of a star-shaped obstacle with a boundary. Assume that a : — t- R 
satisfies for some constants M, i?, oq > 

< a < M, inf a{x) > Oq, q q\ 

f:R^RtsC' satisfying /(O) = f(0) = /"(O) = 0, f{u) = f{\u\), and 

9{u) > 0, < f{u) < Co{\u\^ + g{u)), (1.10) 

for some constant Cq > 0. Then there exists 7 > determined only by N, M, 
R, ttQ and Cq, such that for every initial data {u{0,x),ut{0,x)) with finite energy 
E{u] 0) < 00, there exists a global weak solution u of (11. ip satisfying 

E{u;t) <2e-^'E{u;0) (1.11) 

for all t >0. If (ll.lOp is weakened to 

g{u) > 0, < f{u) < Co(|m|2 + + giu)) (1.12) 

with some 2 < q < 2* := 2N/{N - 2) (2 < q < 00 if N < 2), then ffLTTj) still holds, 
but 7 > depends (non-increasingly) on E{u] 0) and q too. 

Note that the rate 7 is uniform in the whole energy space under the condition 
fll.lOp . while under (11.120 it is uniform on each ball of the energy space. In any 
dimension N > 1, examples of the nonlinearity include all defocusing powers and 
exponentials 

/(m) = Aexp(/i|M|^)|Mp+", (A,/i,i/,a > 0), (1.13) 

and their convex combinations. The global wellposedness in the energy space is so 
far available under the growth condition {H^ subcritical or critical case) 

2N 

If'im) - f\u,)\ < C[\m\ + |«2ir>i -u,\, p+l<2* = (1.14) 

with suitable modifications when 2* < 3 (A^ > 6), and exponential growth condition 
for = 2 (see, e.g. [HI |32l [Ej)- Once we have global wellposedness, the above 
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uniform decay of course applies to the unique solution, but our theorem applies even 
if we do not know unique existence (i.e. in the supercritical case). 

Let us compare our result with the preceding ones for the stabilization problem 
of (11. ip . Zuazua in jlHl SI] first proved the exponential decay by using the unique 
continuation result by Ruiz [31], in bounded and unbounded domains for nonlinear 
term /' mapping into L^, i.e. fll.l4p with p < N/{N — 2), and with a{x) effective 
on the whole boundary and around the infinity. Then Dehman in [3] relaxed the 
nonlinearity condition by making use of the Strichartz estimates, and propagation 
of microlocal defect measures, introduced by Gerard [9J to describe possible lack 
of compactness in bounded sequence of linear solutions. The latter result was 
further extended to the full range of subcritical case, i.e. (11.141) with p + 1 < 2*, 
by Dehman-Lebeau-Zuazua [7]. They used a linearization argument enabling them 
to propagate the possible singularities (lack of compactness) along bi-characteristics 
to the region where the damper is effective. However this linear approximation 
might fail in the critical case p + 1 = 2*, see for example Gerard [TU]. The 
critical case was finally solved by Dehman-Gerard [B] using the nonlinear profile 
decomposition of Bahouri- Gerard [2] instead. 

We point out that these proofs do not explicitly provide any bound on the decay 
rate 7. Recently Tebout [38] gave another proof based on a Carleman estimate 
due to Ruiz pi]. His proof is more quantitative, but also requires the condition 
p < N / {N — 2) to bound the nonlinear term in the Carleman estimate. 

Our proof of Theorem 11.11 gives explicit estimates on 7, and moreover it is very 
elementary, especially compared with the above mentioned arguments. It consists 
mainly of three ingredients: (1) a space-time a priori estimate of the Morawetz type, 

(2) a weighted Sobolev-type inequality adapted to the Morawetz-type estimate, and 

(3) an energy equipartition estimate adjusted to the weight of the Morawetz-type 
estimate. 

The Morawetz space-time integral estimate was originally derived in [21] to have 
some weak decay of local energy for the undamped nonlinear Klein-Gordon equation 



in R^, and then later applied to the nonlinear scattering theory of the undamped 
equation by Strauss- Morawetz [35]. In this paper we use a modified version of the 
estimate derived by Nakanishi [23] for the undamped equation, in the form 




which is available in all spatial dimensions. We emphasize that it is essential to use 
the quadratic term (via the weighted Sobolev inequality), instead of the nonlinear 
term g{u), in order to get decay estimates independent both of the nonlinearity and 
of the energy size. This idea was first used by Nakanishi in [23] for the nonlinear 
scattering theory. 

The energy equipartition property was also used in the preceding works for the 
stabilization, but the mass energy term |mP has prevented them to use it in a 
completely nonlinear way. Introduction of the Morawetz-type weight as well as the 
weighted Sobolev inequality enables us to control that problematic term directly. 



U-AU + U + f\u) = 0, 



(1.15) 





dxdt < CE{u] 0) 



xeneR^, N >l 



(1.16) 
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However, our reliance on the Morawetz-type estimate brings certain limitations of 
our approach. Since it represents the global dispersive property of the solutions, it 
seems difficult to use in bounded domains, or more precisely in presence of trapped 
rays. This is a big disadvantage compared with the preceding works, which equally 
apply to bounded domains (provided that one has enough linear estimates). 

Another, more technical restriction in our theorem is the condition fll.lOp posed 
on g. The preceding works require similar conditions to get the global stabilization 
(namely a decay rate independent of the energy size), but they can also work with 
weaker conditions such as uf'{u) > for local stabilization (decay rate dependent 
of the energy size). 

An example of / which satisfies fll.l2p but not f ll.lOp is given by 

f{u) = \u\Hog{l + \u\^), (1.17) 

i.e. a logarithmic perturbation of a linear term, whereas similar perturbations of 
superlinear powers satisfy fll.lOp . for example 

/'(«) = u^log(l + |nn, (1.18) 

which was considered by Tao in [3^ as an supercritical case in M^, getting global 
solutions for smooth radial initial data (cf. [SUj for the non radial case). 

1.2. Focusing case. In Theorem 1 1.1 1 for the defocusing case, our assumption fll.lOp 
on / is used only to discard nonlinear terms in the Morawetz-type estimates. That 
term can be controlled, if we have the nonlinear scattering for the undamped NLKG 
fll.26p . and a global perturbation argument for the Strichartz norms. The first 
argument provides us with global bounds on the Strichartz space-time norms of the 
undamped equation, while the second one allows us to transfer the global Strichartz 
bounds to the damped equation. The two arguments are indeed satisfied, even for a 
focusing nonlinearity, if we restrict to the solutions with energy less than the ground 
state (a solution of the static equation with possibly a modified mass) and inside 
"the potential well". We should point out, however, that the previous works in 
the defocusing case [HI |5l [TJ Ej could also extend to the same setting, as far as 
the nonlinearity satisfies their growth conditions. See [8] for a result for the wave 
equation on a bounded domain. 

For simplicity of presentation, we consider only the whole space as the domain, 
and typical nonlinearities of one of the following examples, though our argument 
applies to more general nonlinearities. 

(1) > 1 and / is subcritical: We assume 

/(m) = Ai|Mr + --- + Afc|Mr, (1.19) 

for some Ai, . . . , > and 2^ < pi < ■ ■ • < < 2*, where 2* = 2 + A/N is 
the critical power and 2* = 2A^/ (A^ — 2) is the critical one. 

(2) critical case. We assume A^ > 3 and for some A, 

f{u) = \uf. (1.20) 

(3) 2D exponential case. We assume that N = 2 and, 

f(u) = e^"l"l' - 1 - A7r\u\^ - {An\u\y/2. (1.21) 
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Note that the scattering for the lower critical power 2^ is still open for large data 
even in the defocusing case, and the subtracted term in f ll.2ip corresponds to 
that power. 

Global solutions to fll.26p do not always exist in these cases. On the one hand, 
when the data is small in the energy space, it leads to a global solution which 
scatters. On the other hand, when the data is sufficiently large, then the solution 
can blow up in finite time, for example if the energy is negative. The threshold 
between these two dynamics can be constructed as follows. 

Let J be the static energy functional given by 

J{^):= f [\W^\'' + M^]dx-2F{^), F{^):= [ f{^)dx, (1.22) 
and denote its derivative with respect to amplification by 

K{u):=2 \Vu\'^ + \u\'^ -uf'{u)dx. (1.23) 

Consider the constrained minimization problem 

m = mf{J{ip) \ 0^^e H\R'^), K{ip) = 0}, (1.24) 

and define 

/C+ = {{uo,ui) G H\R^) X L\R^) \ E{uo,ui) < m, K{uo) > 0}, 
/C- = {{uo,Ui) G H^R^) X L2(M^) I E{uo,Ui) < m, K{uo) < 0}. 
Then for the undamped NLKG 

uu- Au + u + f'{u) = 0, 

the solutions with energy below m are divided as follows: 

• If {u{0),ut{0)) G /C+, then the solution u of fOHD is global. 

• If {u{0),Ut{0)) G /C^, then the solution u of (11.261) blows up in finite time. 

In the subcritical case, this follows essentially from Payne-Sattinger [28] and the 
local wellposedness in the energy space. 

Recently, Ibrahim-Masmoudi-Nakanishi [IB] extended it to the critical and 2D 
exponential cases. Moreover, they proved that in all the cases (l)-(3), the solutions 
of f ll.26p in /C"*" scatter as t — )■ ±oo with global Strichartz bounds, implementing 
Kenig-Merle's argument [18] to the equation without scaling invariance. In 
addition, they gave the following characterization of m (see [T6l Proposition 1.2]), 
which is again well known in the subcritical case: There exist Q{x) and c G [0, 1] 
such that 

m = r{Q) = J |VQp + c\Q\^dx - 2F{Q), (1.27) 

-AQ + cQ = f'{Q), II Vgili. + cligili. < oo, (1.28) 

and J'^iQ) is the minimum among the solutions of (11.280 . i.e. Q is the ground state 
for (11.280 . c = 1 in the subcritical case f ll.l9p and the exponential case (ll.2ip . and 
c = in the critical case (ll.20p . 

Hence transferring the global Strichartz bounds to the damped NLKG, we obtain 
the following stabilization in the set K^. 



(1.25) 
(1.26) 
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Theorem 1.2. LetN >l,VL = and f{u) given by either (1), (2) or (3). Assume 
that a : — )■ M satisfies (11.91) with some M,R,ao > 0. Let < Eq < m. Then 
there exists 7 > determined by Eq, f , M , R, o-nd N such that for any initial 
data satisfying E{u;0) < Eq and K{u{Q)) > 0, we have a unique global solution u 
of (11.11) with exponential decay 

j^Ep{u- 1) < E{u- 1) < 2e-^'E{u- 0), (1.29) 

for all t > 0. Recall that Ep is the free energy given in (II. 5p . 

In the above Theorem, the dependence of 7 upon Eq is inevitable, since Q is a 
solution of the damped NLKG without any decay. 

In the sequel, A < B 01 B >^ A means A < CB for some positive constant C. 

2. Proof in the defocusing case: Theorem 11.11 

In this section, we prove the uniform exponential decay in the defocusing case, 
i.e. Theorem 11.11 The proof relies solely on the multiplier argument, without any 
estimate for the linearized equation. However in order to justify some energy-type 
computations, we first approximate the nonlinearity so that we have the global well- 
posedness, derive the uniform decay for the solutions of the approximate equations 
with the same initial data, and then take the weak limit to the original equation. 
Note that it is in general hard to get any information for a given weak solution in 
the absence of the wellposedness. 

2.1. Truncation of nonlinearity. First we reduce to the subcritical case by 
approximation of the nonlinearity. Note that this procedure is not entirely trivial. 
On the one hand, we should not simply cutoff the nonlinearity, since g vanishes when 
/' is linear. In other words, approximate nonlinearity must be at least superlinear. 
On the other hand, we need global wellposedness for the approximate equations to 
justify our energy-type computations. The global wellposedness is available in the 
energy subcritical and critical cases, for which growth condition should be imposed 
on /", rather than on / or /'. Consequently, contrary to the classical methods as 
in [331 [ini [31], we need two step approximations of the nonlinearity, one from below 
and another from above. 

Let V{u) = f{u)/\u\'^ be the nonlinear potential energy. Then we have f{u) = 
V{u)\u\'^ and g{u) = V'{u)\u\^u, so our conditions on / are rewritten as 

V(0)=0, uV'{u)>0, V{u) <Co{l + uV'{u)), (2.1) 

and the weaker condition (I1.12p is given by 

V{u) < Co(l + uV'{u) + |m|^-2). (2.2) 

We construct a sequence of approximate nonlinear potentials Vk as follows. First 
we fix ^ G (0, 1), depending only on the dimension A^, such that 2 + < 2*. If there 
exists some A; G N such that for all Izl > 



zV'iz) < kV'{k)\z/kf, 



(2.3) 
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then / is subcritical, i.e. < \zV\z)z\ + < l-zl"*""^^. In this case we 

put fk = /. Otherwise, there exists a sequence of A; G N along which V'{k)k^^^ is 
increasing. For those k we define the first approximation : R — > M by 

Vk{z) = V{z) {\z\ < k), zVliz) = mm{zV'{z), kV'{k)\z/kf) {\z\ > k), (2.4) 

and /fc : M ^ M by 

Mz) = Vk{z)\z\''. (2.5) 

Then clearly we have Vk > 0, > 0, and is locally Lipschitz. Our choice of 
k implies that zVl{z), Vk{z) and fk{z) are all increasing in k at each z. The V 
dominance is inherited by Vk as follows, where the constants may depend also on 6: 
for z > k we have 

Vk{z) < V{k) + mm{ rV\y)dy, f V\k)\y/k\'-'dy) 

Jk Jk (2.6) 

< min(V(z), V{k) + V'{k)\z\'k'-') < 1 + zVliz), 

and similarly for z < —k, while it is trivial for \z\ < k. ( I1.12p is inherited in the 
same way. 

Although the growth for fk and is subcritical, fj! might not be subcritical 
when /" is oscillating for large u. Hence to have global wellposedness for approximate 
equations, we need further approximation V^; : R — t- M defined ioi I > k by 

Vkiiz) = Vk{z) {\z\ < I) zVii{z) = IV^z/lf {\z\ > /), (2.7) 

and /fe; : R — 7- R by fki{z) = Vki{z)\z\'^. Then /^^ is locally Lipschitz and 

\fLi^l)-fMi^2)\<i\z,\ + \z2\Y\z,-Z,\. (2.8) 

Now that the approximate nonlinearity is subcritical (in fact we can make it as 
close to be linear as we want), we have the global wellposedness for 

+ a{x)dtu - Au + u + f'^iiu) = 0. (2.9) 

Now let Uki be the global solution for each k < I with the same initial data, and 
assume that we have 

Ekiiuki; t) < 2e~^'Eki{uki; 0), (2.10) 

for some 7 > 0, all t > 0, and all A; < /, where Eki denotes the energy for the 
approximate equation 

EM{v;t)= I \v\^ + \Vv\^ + \v\^ + fki{v)dx. (2.11) 
Jn 

Then Uki is uniformly bounded in x L^, and so the classical weak compactness ar- 
gument (cf. [SUES]) gives us a subsequence, which we still denote by Uku converging 
as / — !■ 00 to some function Uk{t,x) in the sense 

Ukiit) -> Ukit) in C([0, 00); w-H^{n)), 

iiki ilk in w-L^„^((0, 00); L {Vl)), 

where w-X denotes the space endowed with the weak topology, and also Uki{t^ x) — )■ 
Mfc(t, x) pointwise for almost every a; G at alH > 0. Note that in this limit / — t- 00, 
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the nonlinear energy fki and its limit fk are subcritical, and so are dominated by 
the norm through the Sobolev embedding. Hence the limit Uk is a global weak 
solution to 

d'^u + a{x)dtu -Au + u + f'k{u) = 0, (2.13) 

and moreover, by Fatou's lemma, the uniform decay for Uki, and the Lebesgue 
dominated convergence theorem respectively, we have 

Ek{uk]t) < Wmini Eki{ukut) < liminf 2e"^*Efci(M; 0) = 2e"^*Efc(Mfc; 0). (2.14) 

Next for a subsequence of — )■ oo we have a similar convergence Uk ^ u as above, 
and by Fatou, the uniform decay of Uk and the monotone convergence theorem 
respectively we get 

E{u;t) < \im inf Ek{uk;t) < liminf 2e"^*Efc(M; 0) = 2e~^^E{u]0), (2.15) 

fc— ^-oo fc— >oo ^ ' 

and u is a global weak solution of the original damped NLKG. Note that we cannot 
use the dominated convergence theorem in this step, since the nonlinear energy in 
the limit is no longer controlled by the Sobolev embedding. 

It remains to prove weak convergence of the nonlinear term in the equation. 
Multiplying the equation of Uki with Uki and integrating it in space-time, we get 

[{uki\uki) L^]o + W^UkiWli - \\uki\\l2dt + {fki{uki)\uki)L2dt = 0, (2.16) 
Jo " " Jo 

which implies through the energy identity that 

{fkiiuki)\uki)L2dt<E{u;0). (2.17) 







Then by the dominated convergence theorem for the limit / — ?■ oo, and by the 
monotone convergence theorem for the limit /c — cxd, we get 

/ {fk{uk)\uk)L2dt+ I {f{u)\u)^,dt<E{u;0). (2.18) 
Jo ""Jo 

For any L > 1, let f^{u) = min(/(u), /(L)) and fl^{u) = mm{fkiu)Jk{L)). Then 
for any bounded domain K G Q we have 

\\fk{Uk) - /(u)||li([o,t]xX) < Wftiuk) - /^WILi J[0<t<r]xi^) 

+ \\ukfk{uk)/L\\Li jo<t<T, K|>L) (2.19) 
+ hfiu)/L\\ Ll^{0<t<T, |m|>L)- 

On the right, the first term tends to as A; — )■ oo by the dominated convergence 
theorem, while the other two terms are bounded by E{u;0)/L thanks to (12.181) . 
Hence fkiuk) fiu) in L/q^((0, oo) x VL), and so m is a global weak solution of the 
original equation satisfying the uniform exponential decay. 
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2.2. Uniform decay for subcritical nonlinearity. Now it suffices to prove tlie 
uniform exponential decay in the subcritical (and defocusing) case, with the 
decay rate independent of the nonlinear term and the initial data. Since we have 
global wellposedness in this case, energy-type computations used below are easily 
justified by standard approximation argument. In particular we have the energy 
identity (11. 6p for any finite energy solution u. Denote the energy decrement on the 
interval [0, T] by 

A{u-T):= ! [ aix)\u\^dxdt. (2.20) 



We are going to show that for some quantitative T > and 6 > 0, and for all 
solutions u with finite energy we have 

A{u;T)>6E{u;T). (2.21) 
Then the energy identity implies that 

E{u] t + T) < E{u] t) - 6E{u; t + T), (2.22) 

for any t > 0, and hence by repeated use of it we get the exponential decay (11. lip 
with the rate 

7 = T-Mog(l + 5). (2.23) 

Hence we are mostly concerned with the case where A{u;T) is small compared 
with E{u;0). The main ingredients of our argument are the energy equipartition 
property, and the following Morawetz-type estimate, which is a trivial modification 
of that [21] for the undamped equation. Let 

2 U|2 



(A^-1) t^-\x 



A=|(t,x)| = v/^^TN^, q='-^+ ,3 ' 

2A /2.24) 

—tUt + X ■ vu 

m[u) = h uq. 

A 

Multiplying equation (II. ip with m{u) and integrating it in the light cone, we get 
the Morawetz-type estimate with a space-time weight: 

/ Len A3 ^ ^^^^'^ ^ aum{u)dxdt < E{u] 0) (2.25) 

for any 1 < S <T and any solution u. See [231 El] for the details of the computation 
in the whole space MJ^ . Without loss of generality, we may now choose 

l< R< S = m<T, (2.26) 

where R is the radius for effective damping given in (II. 9p . 

The damping term in (I2.25P is dominated by the energy decrement 

\aum{u)\dxdt < ||ct||//J \/ A{u] T)||m(M) ||^2 (^s<t<T) 

<Mi/^v/I(^Ti/2v/EM) ^^'^^^ 
<E{u-Q) + MTA{u-T). 
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The kinetic part in (12.251) is discarded by 
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T 

s Jn 



xu\ 

"A3" 



-dxdt < 



\x\<R 



<E{u-0) + 



u\^dxdt + 
A{u-T) 

CLqR 



a X m 



dxdt 



(2.28) 



Using that g > as well, we thus get 



\Vu\ 

|x|<t 



+ 



{t^ - \x\^)g{u) 



dxdt < E{u]0)[l + fi{u;T)]. 



where we put 



fi{u; T) := [MT + {aoR)-^]A{u; T)/E{u; 0). 



(2.29) 



(2.30) 



Next, we fix a smooth cut-off function x{^) ^ C^(M^) such that x{^) = 1 for 
\x\ < R and x{^) = for |a;| > 2R, and multiply (II. ip with Integrating it in 
fl, we get 



dt{u\xu) 



Ll 



■ |Vm|2 - \u\^ -ufXu)] dx 

I [auxu + uVu ■ \/x\dx. 
Jn 



(2.31) 



Then by partial integration in t, and using the support property of x, 



s Jn 



\u 



X dxdt 



t 



< 



< 



^ dt{u\xu) 



Ll 



^dt + 



|Vup + |up + \uf'{u)\ + \au\ 



E{u;0) + MA{u]T) + 



xen 

x\<2R 
T 



-dxdt 



(2.32) 



\u\ 



I "/'HI 



xen 

\x\<2R 



t 



dxdt. 



For the last term we use the assumption on / to replace it by 

\uf{u)\ < g{u) + 2f{u) < (1 + Co)giu) + Co\u\\ 



(2.33) 



Then the g term as well as the gradient term is treated by the Morawetz-type 
estimate. For the remaining |mP part, we apply a weighted Sobolev inequality to 
the quadratic part of the Morawetz-type estimate, which implies 



xen t 

'x\<t 



-dxdt < E{u- 0) 



p/2-l 



E{u-Q) + 



\x\<t 



\xu + tVwl 
A3 



-dxdt 



(2.34) 



for all 2 + < p < 2N / [N — 2). The proof was given in for the undamped 
NLKG, which works verbatim for the damped NLKG. Then applying Holder and 
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the support property of x, we get 




\u 



X dxdt < 



s Jn 




\u\ 



xen t 

x\<t 



-dxdt 



2/p 



\x\<2R 



<E{u;0)[l + fi{u;T)f^{R''\ogTy' 



^ dt , 
— dx 

' t 

2/p 



1-2/p 



(2.35) 



On the other hand for |x| > i?, -it is estimated by 




x£Vl t 
\x\>R 



■dxdt < 




a[x]\u 



s Jn 



aot 



-dxdt < ^ ' ^ 



aoS 



(2.36) 



Putting together f l232|) . I^M), (EJSD and I^MI, we obtain 



t 



dt<{l + C^)E{u-Q)[l + ^i{u-T)] 

+ E{u- 0)[1 + ii{u- T)fP{R'' logT)i- 



(2.37) 



2/p 



This is an upper bound. To get a lower bound, we go back to f l2.3ip and replace x 
by 1. Then we obtain 

dt[{u\u)^2 + {au\u)^2/2] = ||?i||^2 - K{u), 



where as before we denote 

K{u) = / \V ul"^ + \u\'^ + uf'{u)dx = E{u;t) - \\u{t) + / g{u)dx 
Jn 

> E{u;t) - \\uit)\\%. 

Hence we get 

2\\u{t)\\l2 > E{u; t) + dt{u\u)^2 + {au\u) j^^, 
and integrating it by dt/t, 



(2.38) 



(2.39) 



(2.40) 



dt 



2\\u 



Ll 



dt 
T 



E{u; t) + dt{u\u)^2 + {au\u)^2^ 

> E{u;T)log- + [{u\u/t)^2fs + / 

s 



MA{u; T) 
T 



^ 2 dt 



(2.41) 



t2 



> E{u; T) log — - E{u; 0)/R - MAiu; T). 



Combining it with f l2.37p . we obtain 

T 



Hence if 



E{u- T) log ^ < (1 + Co)E{u- 0) [1 + ^^{u■ T)] 

+ E{u- 0)[1 + ii{u- T)]2/^'(i?^logT)i-2/P. 

[1 + MT + (aoi?)"^]A(u; T) < E{u; 0)/2, 



(2.42) 



(2.43) 
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then we have T) < 1/2, E{u; T) > E{u; 0)/2 and 

logT < 1 + Co + logi? + iR^\ogTY-''/P, (2.44) 

which imphes 

log T < 1 + Co + log i? + < 1 + Co + . (2.45) 

In other words, there exists an absolute constant C* > such that we have (12.2 ip 
and hence (11.111) with 

logT = C.(l + Co + i?'), 6 = [l + MT+{aoR)-Y^/2, (2.46) 

where we chose p = 2 + 4/A^. 

Finally we consider the weaker assumption (I1.12p . Then (I2.33P is changed to 

\ufiu)\ < (1 + Co)9iu) + Co(|n|2 + {u^. (2.47) 

where without loss of generality, we may assume that q > 2 + A/N. Then to bound 
(I2.32P we need also (12.341) with p = q, and consequently (I2.37P is modified to 

^ y^rft < (1 + Co)Eiu; 0) [1 + /i(«; T)] + C,E{u- ^^^[1 + /i(«; T)] 

+ E{u- 0) [1 + /.(«; T)]2/^(i?^ log T)i-2/^ 
and so the decay rate in this case is given by 

logT = C,(l + Co + i?' + Co^(u; 0)^/2-1), 5 = [I + MT + [aoRy^]-^ /2. (2.49) 
□ 

3. Proof in the focusing case: Theorem 11.21 

In this section, we prove the exponential decay in the focusing case, namely Theo- 
rem [L2l In order to proceed in the same way as in Section I2l2| we need the following 
three ingredients: 

(1) Upper bound on Ep{u; t) in terms of E{u] t). 

(2) Lower bound on K{u(t)) in terms of E{u] t). 

(3) Upper bound on the space-time integral of \g{u)\l\ in terms of the energy 
and Strichartz norms of the undamped equation. 

3.1. Variational properties. The first two are well-known and coming from the 
variational characterization of Q. In fact, since \uf'{u)\ > (2 + 4/A^)|/(-u)|, we have 

(A^ + 2)E{u; t) > NK{u{t)) + 2Ef{u; t) + j N\u\'^dx > 2Ef{u; t), (3.1) 

and for any Eq G (0, J{Q)) there exists z/ > determined by A^, / and Eq such that 
for any G H^{MJ^) satisfying Ji^^) < Eq and K[ip) > we have 

K{^) > lyMli. (3.2) 

For a proof of this lower bound, see [161 Lemma 2.12]. In the exponential case 
/ = fexp {N = 2), we have another important property fl6\ Lemma 2.11] 

\\Vu{t)\\l2 + \\u{t)\\l2 < E{u;t) <Eo<m<l, (3.3) 
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which ensures that the solution stays in the subcritical regime for the Trudinger- 
Moser inequahty 

l|V„|t,.,.., < 1 ^ / - l)dx < . (3.4) 

3.2. Global Strichartz bound from the scattering. The integral bound on 
\g{u)\/X requires transfer of global Strichartz norms from the undamped equation. 
Specifically, we use the following Strichartz norms: 

\MstT ■■= lkllL2*(o,T;<^2(M^)) + II"IIl?+'/''^-''(0,T;B,%(^_^,^{M^))' (3.5) 

if > 3, 

ll'^IUtT ■= W'^'W Lf* {0,T;BI^J^2(^2)) + ll'"llL4(o,r;B;^1,(M2))) (3.6) 

if AT = 2, and 

\H\str ■= \MLf(0,T;Bl(^{R)r (3.7) 

if = 1, where B^^ denotes the inhomogeneous Besov space. 

The main result of [16] asserts that under the assumptions of Theorem 11.21 there 
exists a unique global solution v{t, x) : Mi+^ M for the undamped NLKG (O^ . 
which satisfies 

sup EF{v]t) < A < oo, \\v\\ st^ < B < oo, g\ 

for some A and B determined by iV, / and Eq. Moreover, one has the following 
global perturbation result |T6l Lemma 4.5]: There exists e > 0, determined by N, f 
and Eq, such that for any T > 1 and any u{t, x) : [0, T] x — )■ M satisfying 

Ef{u -V;QY/^+\\u-/\u + U + f\u) |Li(o,T;Li) < ^, (3-9) 

we have 

\\u - v\\stT <C < oo, (3.10) 

where C is also determined by N , f and Eq. 

We apply the above perturbation to the solution u of the damped NLKG (11. ip 
with the same initial data Ep{u — v;0) = 0. Since 

hu\\LUo,T;Li) < M'/'T'/'A{u; Tfl\ (3.11) 

we get 

\\u\\stT <C + B < (X), (3.12) 

provided that 

A{u-T) <{MT)-^e. (3.13) 
Otherwise we get the exponential decay (12.231) with 5 = {MT)~^e / Eq. 
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3.3. Morawetz estimate from the Strichartz bound. By using the Strichartz 
bound (13.121) . we can bound the nonlinear part in the Morawetz estimate 

' Mjf^M, (3.14) 
A 

as foUows. Denote M := [(A^ + l)^o]^^^ + C + B so that we have 

Epiu-tf'^ +\\u\\st^<U. (3.15) 

First we consider the case iV > 3. For the critical power, we get by the Holder, 
Sobolev and interpolation inequalities 

I 1 - III I 1 1 II II 2* 

Jr^ fI ^ ) C3_^gN 

4 2(iV+l) ^ ' ' 

< Il7/I|2* < lU/ll <'^-i'<'^-2' lU/ll < AA^* 

~ ll^llif {0,T;B1/';) ~ II ll^?°(0,T;Hi)ll^ll5tT ^■'^ ' 

where L^''' denotes the Lorentz space on R^. For the critical power, we have 

^ [ w < III I 1 

< lU l|2* < |U ||2* < iV/'2* 

Next we consider the case N = 2 with /p, p > 4 = 2^. Then we have 



T 

I \'l I \^ 

P 



-dxdt < II \x\ ^11 r2,oo IImI 



-/r^ fi 



< II l|P < II l|P"^ II 11^ (3 18) 

~ ll^llLf(0,T;B^/|') ~ ll^llir{0,T;iJi,2)ll^llL^(0,T;Bl^'*) 
~ ll^llLj°=(0,T;//i)ll^ll5tT • 



For the case = 1 with fp and p > 6 = 2^, we have 

— dxdt < ||^^|lioo(o,2..^^)||M||i6(o,T;L6) • (3.19) 

Finally we consider the exponential case / = fexp {N = 2), which requires a bit 
more work. We start with Holder as in the previous cases 

' \9iu)\ 







\X\ 



dxdt<\\\x\ ^11^2,^ II^(m) 11^1(0^^.^2,1), (3.20) 



and to the L^'^ norm we apply the following exponential estimate. Its version 
was proved in 



Lemma 3.1. Let g = uf'{u) — 2f{u) and f{u) = fexp{u) as given in (ll.2ip . Then 
for any A G (0,1), there exist k G (0,1/8), and a continuous increasing function 
Ca '■ [0, oo) — )■ [0, oo), such that for any ip G H^(M?) satisfying 

||V(^|U2 < A (3.21) 

we have, 

lk(^)||L2,i < CA{\Mm)\\v\\%/i-A\v\\%>^/2^ + Il'^|li6,2. (3.22) 



2/k,2 
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Note that f l3.2ip is satisfied by u(t) thanks to the variational property (13.31) . 
and the coefficient CA(||M(t) is uniformly bounded. Thus we can continue the 
estimate (13.201) using the above lemma, 

kMi^™^^ < ||,,||4 II ||4 -L||„||6 

^ 1 1 W 1 1 iO,T;m ) 1 1 W 1 1 ^4 (o,T;Blf2') " " " LUO,T;Bli% ) (3.23) 

It remains to prove the lemma. It is based on a sharp logarithmic inequality 
proved in [131 Theorem 1.3]. 

Lemma 3.2. Let < a < 1. For any A > and any < fi < 1, there exists a 
constant Cx > such that, for any function u G H^{E?) fl C°(M^) 

Ml^ < \\\u\\l^ log (c, + ^>:^Hk^ \ (3^24) 

where := ||Vm||^2 + /^^||w||^2- 

The above logarithmic inequalities constitute a refinement (with respect to the 
best constant) of that appeared in Brezis-Gallouet [1]. 

Proof of Lemma \3. 1\ Without loss of generality, we may assume that A G (1/2, 1). 
Let ip e satisfying (13.211) . We define /i > by 

2fi\l + \\u\\li) = l-A^, (3.25) 

so that we have 

IIv'IIh, <^i - (l + ^')/2< 1. (3.26) 
Now we can choose k G (0, 1/8) and A > 2/7r, depending only on Ai, such that 

27rA(l/4 -k)>1 = 27rA(l/4 + k/2)AI. (3.27) 
For instance, we can choose k = (1 — Al)/8. Note that 

<I^IV^'"''-1)- (3.28) 
In particular, if < Ai, then \g{(p)\ < Iv^l^ and thanks to Sobolev embedding 

\\9{^)\\l^-' ^ Il'^llii2,6 < ||'^||ii2,2 < Ilv^ll|v6- (3.29) 
If ||<y3||Loc > Ai, then we have 

lU^^MI ^ ll^47r|(p|2 II l/2-K|| 47r|mp n 1/2+k n n 4 

||^(<^)||l2,i < ||e I'^i -l|Li ||e -1|Il- ||V5|Il4/.,4/(i/2+k) 

< _ 1 ||l/2-«|| 47r|^P .||l/2+«|| ||4 

— Il'' IIli II -'^IIl°° llrllL4/K,2- 
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The first factor on tlie riglit liand side in tlie above inequality is bounded by tlie 
sliarp Trudinger-Moser inequality. Since ||v9||c'i/4-k > A > 1/2 and Cx > 1, then the 
second factor is bounded by using the logarithmic inequality f l3.24p 



27r(l+2re)A||93|p 

|g47r|¥>|2 _ ]^||l/2+K ^ g27r(l+2K)j|(p||2^ ^ 



(3.31) 

Thus, using Sobolev embedding we obtain 

II^MIIl^.i < ||v?||^i/6 + C7^(||</'||//i)ll^llci/4-.||v^ll^./2 • (3.32) 



2/k,2 



□ 



3.4. Concluding the proof of Theorem 11.21 Thus we have so far either 

A{u;T) > {MT)-^e, (3.33) 

or 

r [ ^-^^dxdt < C{M)E{u- 0) < oo, (3.34) 
Ji Jm ^ 

where C{Af) > is a certain uniform bound in terms of Eq, f and N. In the latter 
case, we obtain from f l2.25p in the same way as in Section ?I72\ using the free energy 
bound (13. ip as well, 

f |VuP + |a;M + tVuP \g(u)\ , , 

/ ^ 4 + -^^\^dxdt 

s J\x\<t ^ ^ (3.35) 

<Eiu;0)[l + ^xiu;T)+Ci^^)], 

ioT 1 < R < S = 3R < T, where fi is as defined in (I2.30p . Then proceeding in the 
same way as before, we obtain 

y^dt < (1 + Co)E{u; 0) [1 + fi{u; T) + C{U)] 

+ E{u- 0)[1 + ii{u- T) + C{N)f^{R'' logT)i-2/f , 

instead of fl2.37p . Further tracking the same way, we replace fl2.39p by the following 
bound due to (13. 2p 

K{u) > u\\umm > iy[E{u-t) - \Ht)\\%], (3.37) 

hence we have 

(1 + iy)\\u{t)\\l2 > uE{u;t) + dt{u\u)L2 + (aM|M)^2. (3.38) 
Integrate it by dt/t, we obtain the following instead of (I2.4ip : 

/■^ df T 

-(1 + v) \\u{t) \\l2 > vE{u- T) log — - E{u- 0)/R - MA{u- T). (3.39) 
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Combining it with f l3.36p . we obtain 

vE{u- T) log ^ < (1 + C,)E{u- 0) [1 + ^^{u■ T) + C{M)] 

+ E{u- 0) [1 + ii{u- T) + C{U)fl^{R'' log Tf-'^'K 

Hence if 

[1 + MT+ {aoR)-^]A{u] T) < E{u; 0) /2, (3.41) 
then we have /i < 1/2, E{u; T) > E{u; 0)/2 and 

u\ogT< log/? + (1 + Co)[l + CiAf)] + [1 + C(Ar)]2/f(i?^ logr)i-2/p^ (3_42) 
which implies that 

logT < + C{Af))il + Co + R'), (3.43) 

where we have chosen p = 2^ = 2 + 4/A^. Thus in conclusion there exists an absolute 
constant C* > such that we have (11.111) with 

logT = + C{^)){1 + Co + R'), 

1 1 1 (3.44) 

S = mm{[l + MT+{aoR)-T'/2,{MT)-h/Eo), 

where the last factor takes care of the case (13.331) . □ 
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